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ABSTRACT 

We compute the luminosity function (LF) and the formation rate of long gamma ray bursts (GRBs) 
by fitting the observed differential peak flux distribution obtained by the BATSE satellite in three 
different scenarios: i) GRBs follow the cosmic star formation and their LF is constant in time; ii) 
GRBs follow the cosmic star formation but the LF varies with redshift; iii) GRBs form preferentially 
in low-metallicity environments. We find that the differential peak flux number counts obtained by 
BATSE and by Swift can be reproduced using the same LF and GRB formation rate, indicating 
that the two satellites are observing the same GRB population. We then check the resulting redshift 
distributions in the light of Swift 2-year data, focusing in particular on the relatively large sample 
of GRBs detected at z > 2.5. We show that models in which GRBs trace the cosmic star formation 
and are described by a constant LF are ruled out by the number of high-z Swift detections. This 
conclusion does not depend on the redshift distribution of bursts that lack of optical identification, 
nor on the existence of a decline in star formation rate at z > 2, nor on the adopted faint-end of the 
GRB LF. Swift observations can be explained by assuming that the LF varies with redshift and/or 
that GRB formation is limited to low-metallicity environments. 
Subject headings: gamma-ray: burst - stars: formation - cosmology: observations. 



1. INTRODUCTION 

Gamma Ray Bursts (GRBs) are powerful flashes of 
high-energy photons occurring at an average rate of a few 
per day throughout the universe. Even though they are 
highly transient events very hard to localize, they are so 
bright that they can be detected up to very high redshift 
(the current record is z = 6.29). The energy source of a 
GRB is believed to be associated to the collapse of the 
core of a massive star in the case of long-duration GRBs, 
and due to merger- or accretion-induced collapse for the 
short-hard class of GRBs (see Meszaros 2006 for a recent 
review). In this paper, we limit our analysis to the class 
of long-duration GRBs. 

One of the main goals of the Swift satellite (Gehrels et 
al. 2004) is to trackle the key issue of the GRB lumi- 
nosity function (LF). Unfortunately, although the num- 
ber of GRBs with good redshift determination has been 
largely increased by Swift, the sample is still too poor 
(and bias dominated) to allow a direct measurement of 
the LF. Many studies (e.g. Lamb & Reichart 2000; Por- 
ciani & Madau 2001 (PM0I); Schmidt 2001, Choudhury 
& Srianand 2002; Firmani et al. 2004; Guetta, Piran 
& Waxman 2005; Natarajan et al.2005; Daigne, Rossi & 
Mochkovitch 2006) tried to constrain the GRB LF under 
the assumption that GRBs trace the observed star forma- 
tion rate, as suggested by the association of long GRBs 
to the death of massive stars. Following these works and 
assuming the most recent star formation rate determina- 
tion, we derive the LF and formation rate of GRBs by 
fitting the observed BATSE differential peak flux distri- 
bution in three different scenarios: i) GRBs follow the 
cosmic star formation and have a constant LF; ii) the 
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GRB LF varies with redshift; iii) GRBs form in low- 
metallicity environments. We check the results against 
the 2-year Swift data, focusing in particular on the large 
sample of high redshift (z > 2.5) GRBs detected by this 
instrument. 
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2. BASIC EQUATIONS 
The observed photon flux, P, in the 



< E < E n 



energy band 
emitted by an isotropically radiat- 



ing source at redshift 
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where S(E) is the differential rest-frame photon lumi- 
nosity of the source, and di,(z) is the luminosity dis- 
tance. To describe the typical burst spectrum we adopt 
the functional form proposed by Band et al. (1993), i.e. 
a broken power-law with a low-energy spectral index a, 
a high-energy spectral index (3, and a break energy Et- 
In this work, we take a = — 1 and (3 = —2.25 (Preece 
et al. 2000), and E b = 511 keV (PM0I). Moreover, 
it is customary to define an isotropic equivalent intrin- 
sic burst luminosity in the energy band 30-2000 keV as 
L = f™°y V ES(E)dE. Given a normalized GRB LF, 
4>(L), and the detector efficiency, e(P), the observed rate 
of bursts with peak flux between P\ and P2 is 
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where dV(z)/dz = ATrcd 2 L (z)/[H(z)(l + z) 2 ] is the comov- 
ing volume element 3 , and H(z) = Hq[Q,m (1 + z ) 3 + ^A + 

3 We adopted the 'concordance' model values for the cosmolog- 
ical parameters: h = 0.7, Q m = 0.3, and Qa = 0.7. 
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(1 — f2jvr — ^a)(1 + ^) 2 ] 1 ^ 2 - Af2 s is the solid angle covered 
on the sky by the survey, and the factor (1 + z)~ x ac- 
counts for cosmological time dilation. Finally, ^gkb{z) 
is the comoving burst formation rate. In this work, we 
assume that the GRB LF is described by 
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3. MODELS 



We consider three different scenarios. In the first 
one, the GRB formation rate is proportional to the cos- 
mic star formation rate (SFR), ^(z), i.e. ^grb(-z) = 
fcGRB^f '*{z) , and the LF does not evolve with redshift, i.e. 
L C \it — const = Lq. The factor kQRB gives the number 
of GRBs formed per solar mass in stars and has units 
of Mq 1 . (in units of M© Mpc~ 3 yr" 1 ) is com- 

monly parameterized with the form proposed by Cole et 
al. (2001) as 



(ai + a 2 z)h 
1 + (z/a 3 ) a * 



(4) 



Recently, Hopkins & Beacom (2006) have provided the 
values of the coefficients a by fitting the available UV 
and far-infrared measurements for z < 6, corrected for 
dust obscuration. In this paper, we adopt their best fit 
parameters: ai = 0.017, a 2 = 0.13, and a 3 = 3.3 (Hop- 
kins & Beacom 2006). The value of 0,4 = 4.3 is taken to 
be slightly lower than the original one in order to match 
the decline of the SFR with (1 + z)~ 3 - 3 at z > 5 sug- 
gested by recent deep-field data (see Stark et al. 2006 
and references therein). 

In the second scenario, while the GRB formation rate 
is still proportional to the observed SFR, the cut-off 
luminosity in the GRB LF increases with redshift as 
L C ut = Lq(1 + z) s . Lloyd-Ronning, Fryer & Ramirez- 
Ruiz (2002), using GRB redshifts and luminosities de- 
rived from the luminosity- variability relationship, found 
that the data imply 5 ~ 1.4 ± 0.5, and we adopt this as 
fiducial value. 

Finally, we consider a case in which GRBs form only 
in environments with metallicity below a given thresh- 
old (no evolution in the LF is considered). In fact, 
some theoretical models (see Meszaros 2006 and refer- 
ence therein) require that GRB progenitors should have 
metallicity < 0.1 Z Q . Observations of GRB host galaxies 
(see Savaglio 2006 and reference therein) seems in agree- 
ment with this prescription, showing that GRB prefer- 
entially originates in low-metallicity regions. Langer & 
Norman (2006) have quantified the amount of star for- 
mation at a given metallicity, using a recent determina- 
tion of the stellar mass function (Panter et al. 2004) 
and the observed mass-metallicity correlation (Savaglio 
et al. 2005). Adopting the metallicity redshift evolution 
derived from emission line studies (Kewley & Kobulnicky 
2005), the fractional mass density belonging to metallic- 
ity below a given threshold, Zth, can be computed as 
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where F (r) are the incomplete (complete) gamma func- 
tion, and T(0.84) ~ 1.122. The GRB formation rate is 



then given by *grb(^) = k G B.B^(z)^*(z). The main 
effect of this convolution is that the GRB formation rate 
peaks at higher redshift with respect to the cosmic SFR. 
We adopt Z t h = 0.1 Z as fiducial value, and, in this 
case, the GRB formation peaks at z ~ 3.5. 

4. GRB NUMBER COUNTS 

The free parameters in our model are the GRB forma- 
tion efficiency ^grb, the cut-off luminosity at z = 0, Lq, 
and the power index, £, of the GRB LF function. Fol- 
lowing PM01, we optimized the value of these parameters 
by x 2 minimization over the observed differential number 
counts in the 50-300 keV band of BATSE. We use the 
off-line BATSE sample of Kommers et al. (2000), which 
includes 1998 archival ("triggered" plus "non-triggered") 
bursts, and for which the detector efficiency is well de- 
scribed by the function e(P) = 0.5[(1 + erf(-4.801 + 
29.868P)] (Kommers et al. 2000). We report the best- 
fit parameters for our fiducial models in Table [1] In the 
last column, we give the reduced x 2 f° r the best-fitting 
model, showing that it is always possible to find a good 
agreement with the data 4 . Note that for the metallicity 
evolution scenario a higher GRB formation efficiency is 
required, since GRBs form only in a (small) fraction of 
star forming galaxies. 

We can now use the best-fit parameters to compute the 
expected differential peak flux distribution of GRBs in 
the 15-150 keV band of the Burst Alert Telescope (BAT) 
instrument onboard of Swift. The results are plotted in 
Figure [1] and compared with the observed Swift/BAT 
data points. All models show a good agreement with the 
data without the need of any change of the GRB LF and 
formation efficiency, indicating that BATSE and Swift 
are observing essentially the same population of GRBs. 
This conclusion is rather insensitive to 20% variations of 
the adopted GRB spectrum parameters, i.e. for the large 
majority of burst spectra (Kaneko et al. 2006). 

5. GRB REDSHIFT DISTRIBUTION 

Our model allows us to compute the expected red- 
shift distribution of GRBs detected by Swift. We de- 
cide to avoid the comparison between model results and 
the overall observed distribution of bursts with known 
redshift, since this procedure implicitly assumes that the 
observed sample of GRBs with redshift determination is 
representative of all detected sources. Moreover, impor- 
tant information are missed by this kind of analysis: for 
example, that many bright GRBs are identified at high 
redshift. So, we try to answer this simple question: is the 
redshift distribution consistent with the number of Swift 
detections at z > 2.5 and z > 3.5? 

The cumulative number of GRBs, identified during the 
two years of the Swift mission at z > 2.5 (left panel) and 
z > 3.5 (right panel), is plotted in Figure [2] toghether 
with model predictions. Note that Swift detections are 
to be considered as a strong lower limit, since many high- 
z bursts can be missed by optical follow-up searches. 

The model with no LF evolution clearly underesti- 
mates the number of high redshift GRB detections at 
any photon flux and no bright GRBs are predicted for 
z > 3.5. We checked that variations of the shape of the 

4 Note that strong covariance on Lo and £ is observed in the pa- 
rameter space surrounding the best-fit parameters (see also PM01) 
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SFR do not affect this result: even assuming a constant 
SFR at z > 2, the model predictions do not change sig- 
nificantly. In fact, for relatively bright GRBs, the rapid 
decline in the LF strongly hampers the detection of high 
rcdshift bursts. Furthermore, our analysis does not de- 
pend on the faint-end of the LF: increasing the popula- 
tion of faint GRBs would decrease the number of high-z 
detections, strengthening our conclusion. So, models in 
which GRBs trace the cosmic SFR and are described by a 
constant LF, are ruled out by the large sample of high-z 
Swift GRBs. 

The number of high-z Swift identifications can be jus- 
tified assuming that the LF varies with redshift. In 
this case, high-z GRBs are typically brighter than low- 
z ones, so that are much easely detected. Assuming 
that the luminosity increases as (1 + z) 1,4 , we find many 
sources at z > 2.5, but the model is barely consistent 
with the number of bright GRBs at z > 3.5. Since 
some high-z sources can be missed by optical follow-up 
searches, an even stronger evolution might be required 
to explain the data. 

Finally, we consider the possibility that GRB forma- 
tion is restricted to low-metallicity environments. In this 
case, the peak of the GRB formation is shifted towards 
higher redshift, so that the probability of high-z detec- 
tions increases. Assuming Z t h = 0.1 Zq, Swift identifi- 
cation are exceeded both at z > 2.5 and z > 3.5 without 
requiring any evolution in the LF. Thus, the model is 
consistent with a fraction of high redshift bursts missed 
by optical follow-up searches. Increasing the threshold 
metallicity will decrease the number of sources at high-z: 
for Zth ~ 0.4 Zq the model becomes inconsistent with 
the number of observed GRBs at z > 3.5. Higher thresh- 
old values would require evolution of the GRB luminosity 
and/or a more gentle decline of the SFR at high redshift. 

In conclusion, the existence of a large sample of bursts 
at z > 2.5 in the Swift 2-year data imply that GRBs have 
experienced some kind of evolution, being more luminous 
or more common in the past. 

6. GRB RATE AT REDSHIFT LARGER THAN SIX 

The discovery of GRB050904 at z = 6.29 (Antonelli 
et al. 2005; Tagliaferri et al. 2005; Kawai et al. 2006) 
during the first year of the Swift mission has strengthened 
the idea that many GRBs should be observed out to very 
high redshift (e.g. Natarajan et al. 2005; Bromm & 
Loeb 2006; Daigne et al. 2006). Unfortunately, no other 
source at z > 6 has been detected in the second year of 
observations. 

In Figure [3l we plot the Swift detection rate expected 
for the three scenarios here considered. Models without 
evolution predict almost no sources to be detected at very 
high redshift. If luminosity evolution (S = 1.4) is allowed, 
~ 2 bursts/yr should lie above z ~ 6 for P > 0.2 ph 
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cm~ 2 s _1 , whereas, in the metallicity evolution scenario 
(Z t h — 0.1 Zq), we expect ~ 8 GRBs/yr, one or two 
being at z > 8. 

The detection rate are found to decrease rapidly with 
increasing peak fluxes. Indeed, it is interesting to note 
that GRB050904 was relatively bright, being its observed 
photon flux P = 0.658 ph cm~ 2 s _1 . At this limit, only 
~ 1 (2) bursts/yr would be at z > 6, if luminosity (metal- 
licity) evolution is assumed. Thus, the lack of very high 
redshift identification in the 2nd year of the Swift mission 
might be due to practical difficulties in the optical follow- 
up of faint GRBs. In fact, no GRB with observed photon 
fluxes below 0.5 ph cm' 2 s _1 has UVOT detection and 
only in a couple of cases a reliable redshift determination 
was possible. So, the identification of just one burst at 
z > 6 in two years of Swift mission is not very surprising. 
On the contrary, the discovery of GRB050904 may sug- 
gest that the Swift follow-up procedure is working very 
well, at least for relatively bright bursts. 

7. CONCLUSIONS 

We have computed the luminosity function and the for- 
mation rate of long GRBs by fitting the BATSE differ- 
ential peak flux number counts in three different scenar- 
ios: i) GRBs follow the cosmic star formation and have 
a redshift-independent LF; ii) the GRB LF varies with 
rcdshift; iii) GRBs are associated with star formation in 
low-metallicity environments. In all cases, it is possible 
to obtain a good fit to the data by adjusting the model 
free parameters. Moreover, using the same LF and for- 
mation rate, it is possible to reproduce both BATSE and 
Swift differential counts, showing that the two satellites 
are observing the same GRB population. 

We have then computed the expected burst redshift 
distribution, testing the results against the number of 
high redshift GRBs, detected during the two years of the 
Swift mission. We find that models where GRBs trace 
the SFR and are described by a constant LF largely un- 
derestimate the number high-z GRBs detected by Swift. 
This conclusion does not depend on the redshift distri- 
bution of burst lacking of optical identification, nor on 
the existence of a decline in the SFR at z > 2, nor on 
the adopted faint-end of the LF. Alternatively, we find 
that the observed number of high-z detection can be 
justified by assuming that the GRB luminosity increases 
with redshift and/or that GRBs preferentially form in 
low-metallicity environments. 

Finally, we have estimated the detection rate of bursts 
at very high redshift. We find that - 2 (8) GRBs/yr 
should be observed at z > 6, if luminosity (metallicity) 
evolution is assumed. The majority of these sources is 
faint and may be missed in optical follow-up searches, 
but ~ 1 (3) GRB/yr should be relatively bright, with an 
observed photon flux in excees to 0.5 ph cm -2 s _1 . 
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TABLE 1 

Best fit parameters for different models. Errors are at la level. 



Model fccWao^M- 1 ) L /(10 51 ergs" 1 ) £ 



no evolution 1.14 ±0.07 9.54 ±4.55 3.54 ± 0.78 0.83 

luminosity evolution (5 = 1.4) 1.05 ± 0.05 0.77 ±0.13 2.19 ±0.95 0.80 

metallicity evolution (Z th = 0.1 Z Q ) 10.0 ± 0.5 16.7 ± 5.7 2.94 ± 0.34 0.84 
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Fig. 1. — DiflFerential number counts for Swift in the 15-150 keV band as a function of the observed photon flux P. The points show the 
observed counts and their Poisson uncertainties (horizontal error bars denote bin size). Dotted lines refers to the model without evolution, 
short-dashed line to the luminosity evolution model (<5 = 1.4), and long-dashed line to the model with the metallicity threshold for GRB 
formation (Z^ = 0.1 Zq). A field of view of 1.4 sr for Swift/HAT is adopted. 
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Fig. 2. — Cumulative number of high redshift GRBs at z > 2.5 (right panel) and at z > 3.5 (left panel) as a function of the observed 
photon flux P in the 15-150 keV band. The number of sources detected in the two years of Swift mission is shown as solid histogram, 
whereas model results are shown with lines as in the previous figure. Note that the observed detections are lower limits, since many high-z 
GRBs can be missed by optical follow— up searches. A field of view of 1.4 sr for Swift/HAT is adopted. 
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Fig. 3. — Cumulative rate of z > 6 GRBs detectable by Swift as a function of the photon flux P. A field of view of 1.4 sr for SwiftfBAT 
is adopted. Lines as in Figure [T] 



